Dimensional reduction of a binary Bose-Einstein condensate in mixed dimensions 
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We present effective reduced equations for the study of a binary Bose-Einstein condensate (BEC), 
where the confining potentials of the two BEC components have distinct asymmetry so that the 
components belong to different space dimensions as in a recent experiment [G. Lamporesi et at, Phys. 
Rev. Lett. 104, 153202 (2010)]. Starting from a binary three-dimensional (3D) Gross-Pitaevskii 
equation (GPE) and using a Lagrangian variational approach we derive a binary effective nonlinear 
Schrodinger equation with components in different reduced dimensions, e. g., the first component 
in one dimension and the second in two dimensions as appropriate to represent a cigar-shaped 
BEC coupled to a disk-shaped BEC. We demonstrate that the effective reduced binary equation, 
which depend on the geometry of the system, is quite reliable when compared with the binary 3D 
GPE and can be efficiently used to perform numerical simulation and analytical calculation for the 
investigation of static and dynamic properties of a binary BEC in mixed dimensions. 

PACS numbers: 03.75.Mn, 03.75.Hh, 3.75. Kk 



I. INTRODUCTION 



It is now well established that, at ultralow tempera- 
ture, dilute Bose-Einstein condensates (BECs) can be ac- 
curately described by the three-dimensional (3D) Gross- 
Pitaevskii equation (GPE) [HQ. Some years ago it has 
been found that, starting from the 3D GPE and using 
a Gaussian variational approach H|, it is possible to 
derive effective one-dimensional (ID) or two-dimensional 
(2D) wave equations which describes the dynamics of 
quasi-lD or quasi-2D BECs as appropriate for cigar or 
disk shapes 0. For instance, in that derivation of the 
effective ID equation the trapping potential is harmonic 
and isotropic in the transverse direction and generic in 
the axial one. The effective ID equation, which is a time- 
dependent nonpolynomial Schrodinger equation (NPSE) 
[f| , has been used to model a cigar-shaped condensate by 
many experimental and theoretical groups [q-tioj. Sim- 
ilar reduction scheme also exists for a disk-shaped BEC 
leading to a 2D NPSE. (The lowest-order approximation 
of the NPSE leads to a GP-typc nonlinear Schrodinger 
equation (NLSE) equation with cubic nonlinearity. Dif- 
ferent variants of this reduction scheme have later been 
suggested, which offer certain advantages in some cases 
Moreover, by using the ID NPSE, analytical and 
numerical solutions have been found for solitons and vor- 
tices [ll[. Recently, a discrete version of the ID NPSE 
(ID DNPSE) has been obtained The ID DNPSE 

has been used to model a BEC confined in a combina- 
tion of a cigar-shaped trap and a deep optical lattice 
acting in the axial direction More recently the 2D 
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NPSE, obtained for a disk-shaped BEC subject to strong 
confinement along the axial direction, has been used to 
predict the density profiles and dynamical stability of 
repulsive and attractive BECs with zero and finite topo- 
logical charge in various planar trapping configurations, 
including the axisymmetric harmonic confinement and 
ID periodic potential [l3||. The reduction scheme has 
also successfully been applied 01 EH to the study of a 
density-functional formulation |15| for a Fermi supcrfluid 
at unitarity as well as to the study of Anderson localiza- 
tion in a BEC 

In a very recent experiment Lamporesi et al. fl7j inves- 
tigated the mixed-dimensional scattering occurring when 
the collisional atoms of two species of a binary BEC live 
in different space dimensions. This achievement opens 
the way to the experimental and theoretical studies of 
a binary BEC with the components belonging to mixed 
dimensions. In the present paper we consider the di- 
mensional reduction of a binary 3D GPE appropriate for 
the study of a binary BEC in mixed dimensions. Starting 
from a binary 3D GPE we derive a binary nonpolynomial 
Schrodinger equation (BNPSE) in reduced dimensions by 
using a Lagrangian variational approach. (The BNPSE 
represent a set of coupled nonlinear differential equations 
with nonpolynomial non-power-law nonlinearities.) We 
also present a simplified version of the reduction scheme 
first, where the reduced equations in the form of a binary 
NLSE (BNLSE) have simpler power-law nonlinearities. 
The BNPSEs depend on the geometry of the system and 
also on the inter-species interaction. In fact, the main 
difference with respect to the case of a single BEC is the 
coupling between the two BECs, which gives rise, quite 
generally, to a integral term in the reduced wave equa- 
tions. (In the past, a binary ID NPSE has been derived 
and used [l|[ to study a binary cigar-shaped BEC, rather 
than a binary BEC in mixed dimensions as considered 
here.) 
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We use the BNLSE and BNPSE to study numerically 
the static and dynamic properties of a 1D-2D binary BEC 
where the first quasi-lD component lies along the z axis 
and the quasi-2D component in the (x, y) plane. For 
small nonlinearities both the BNLSE and BNPSE yielded 
results for density in good agreement with the full 3D 
result. At larger values of nonlinearities the BNPSE pro- 
vided better approximation to the full 3D result. We also 
studied small oscillations of the 1D-2D binary BEC initi- 
ated by a sudden change of the radial angular frequency 
of the quasi-2D component. The BNPSE provided a bet- 
ter description of this oscillation over the BNLSE, when 
compared with the full 3D model, specially for large non- 
linearities and large times. 

The reduced effective BNLSEs with simple power-law 
nonlinearity, but with different dimensionalities for the 
components, are presented in Sec. |TIJ These equations 
are derived by considering the ansatz for the wave func- 
tion that sets the transverse spatial dependence as the 
Gaussian ground states in transverse harmonic traps. 
The transverse spatial dependence is finally integrated 
out and the reduced equations derived in the usual fash- 
ion using a Gaussian variational approach. In Sec. IIIII 
the BNPSEs with different dimensionalities for the com- 
ponents are derived by taking the ansatz for the wave 
function that sets the transverse spatial dependence as 
generic Gaussian states which are not the ground states 
in transverse harmonic traps. The reduced effective 
BNPSEs arc then derived by a Lagrangian variational ap- 
proximation which determines the widths of the Gaussian 
functions in a self-consistent fashion. In Sec. IIV Al we 
present the numerical results, where we compare the den- 
sities and chemical potentials of the 3D binary GPE with 
those of the BNLSE and BNPSE. The binary BEC has a 
disk-shaped component along the z direction coupled to 
a cigar-shaped component in the (x, y) plane. For smaller 
nonlinearities of the 3D binary GPE, we establish good 
agreement between the results for density and chemical 
potential of the effective BNLSE and the full 3D binary 
GPE. For larger nonlinearities better result was obtained 
with the BNPSE. We also studied small oscillations of 
the binary 1D-2D BEC initiated by a sudden change of 
the radial angular frequency of the quasi-2D disk-shaped 
BEC using the BNPSE, BNLSE and the full 3D equa- 
tions. Again the results obtained with the BNPSE were 
in better agreement with those of the BNLSE. Finally, a 
summary of our investigation is presented in Sec [V] 



II. THEORETICAL FORMULATION FOR THE 
BNLSE 



will be assumed to have a quasi-lD or quasi-2D config- 
uration. The quasi-lD shape emerges when the traps 
in the two transverse directions are much stronger than 
that in the linear direction. The quasi-2D shape emerges 
when the trap in the direction perpendicular to the 2D 
plane is much stronger than those in the 2D plane. In 
this fashion one could have an effective 1D-3D, and 2D- 
3D mixed configuration for the traps acting on the binary 
BEC described by effective 1D-3D and 2D-3D BNLSEs. 
There remains one distinct case where one of the compo- 
nents have a quasi-lD shape and the other component a 
quasi-2D shape due to specific trap symmetry of the two 
components. In that case the binary mixture is described 
by a effective 1D-2D BNLSE. We discuss these possibili- 
ties in the following and illustrate the corresponding re- 
duction schemes. 

The action of a system of two BECs is given by 
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where the Lagrangian density is 
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with Uj(r) the external potential acting on the j-th 
BEC component (j = 1,2), which is described by the 
macroscopic wave function i(^j(r, t). Here dt denotes time 
derivative, gj = Airh aj/rrij, g\2 = 2irh 2 ai2/mR, and 
f \ifjj(r,t)\ 2 — Nj, where ctj is the intraspecies scatter- 
ing length for species j and an is the interspecies scat- 
tering length, rrij is the mass of species j, and trr = 
rri\mil (m\ + m-i) is the reduced mass. 

By extremizing the action A with respect to ij>*(r,t) 
and ■j/'I (r, t) we obtain the following binary GPE [19[ 
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These two coupled partial differential equations are both 
in 3+1 (space-time) dimensions. In the general case, 
these can be solved numerically in a routine fashion 
[20I |2lj . However, in many situations of interest, when 
one or both of the BEC components have extreme dis- 
tinct cigar and disk shapes and when we are interested in 
density and dynamics in axial and radial directions, re- 
spectively, the following reduction procedures are of par- 
ticular use. 



We shall now explicitly demonstrate the reduction pro- 
cedure in three different types of binary mixtures: 1D- 
3D, 1D-2D, and 2D-3D. We shall consider a binary BEC 
where the trap anisotropy of the two components are 
distinct. When the component 2 has a fully asymmetric 
trap, the component 1, due to its specific trap symmetry, 



A. 1D-3D binary BEC 

Let us suppose that the confining potential of the first 
BEC is 

EMr) = \m^\ x x 2 + u 2 y y 2 ) + V^z), (5) 
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where V\(z) is a generic potential in z direction, and 
uj\ x and uiy are the angular frequencies of the harmonic 
trap along transverse x and y directions. The transverse 
harmonic potentials are assumed to be so strong that the 
first BEC is quasi-lD, i.e. in an elongated cigar-shaped 
configuration. 

We can then use the variational wave function [fl 



-x 2 /(2a 2 lx ) 



-V 2 /(2al y ) 



for the first BEC, where a-y x = yjh/ (m\u>i x ) and a± y = 
yjh/ {miLOiy) are the characteristic harmonic lengths of 
the oscillators in x and y directions, respectively. Conse- 
quently, the action of the system becomes 



effective BNLSE for fi(z,t) and tp2(p,z,t) 
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where a\ p = a\ x = a\ y = y/h/Jm^Uhp). In the coupled 
set of Eqs. (fl"3|) and (| 14[) . the first is intcgro-diffcrcntial 
in 1+1 dimensions and the second is differential in 3+1 
dimensions (but in 2+1 variables). 



A = \ L\ dzdt 
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We now extremize the action A with respect to f*(z,t) 
and $j[(r,i), to obtain the effective BNLSE for fi(z,t) 
and ip2(r,t): 
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Equations (|lip and (fT2")) are simplified if we impose 
cylindrical symmetry in the harmonic potential acting 
on the first BEC, i.e. u>\ x = u>\ y = U)\ p , and also 
in the generic potential acting on the second BEC, i.e. 
[/ 2 (r) = U 2 {p,z) with p = (x 2 + y 2 ) 1 ! 2 . In this way we 
can choose ^> 2 (r,i) = ^ 2 (p, z,i) and find the following 



B. 1D-2D binary BEC 

Parallel BECs: Let us suppose that the confining 
potential of the first BEC is still given by Eq. ([S]), while 
the confining potential of the second BEC is 



&2 (r) = ^m 2 wl x x 2 + V 2 (y,z) 



(15) 



where V2(y,z) is a generic potential acting in the (y,z) 
plane. In addition we impose that the harmonic po- 
tential of angular frequency UJ2 X along the x axis is so 
strong that the second BEC is quasi-2D, i.e. a very flat 
two-dimensional configuration in the (y, z) plane. Con- 
sequently, we have binary 1D-2D BEC where the first 
BEC is cigar-shaped along the z axis and the second disk- 
shaped in the (y, z) plane. 

We can then use the variational wave function of Eq. 
© for the first BEC, while for the second BEC we adopt 
the following variational wave function Q 



e -x 2 /(2a 2 2 J 
^(r.t) = -i—r^ rrrpr <f) 2 {y,z,t) , 
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where a 2x = y/ ft/(m 2 w 2a; ) is the characteristic harmonic 
length of the oscillator in x direction acting on the second 
BEC. In this way the action of the system becomes 



A — J L\ dzdt + / £ 2 dydzdt + / £ 12 dydzdt , (17) 
where L\ is still given by Eq. ((8j), while 
£2 
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We now extremize the action A and find the following 
effective BNLSE for fi{z,t) and (/> 2 (y, z,t); 
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These arc two coupled equations, where the first is 
intcgro-differential in 1+1 dimensions and the second is 
differential in 2+1 dimensions. Equations (|2 1 [) and (|22|) 
describe the binary system where the axis of the cigar- 
shaped BEC along z direction is in the (y, z) plane of the 
two-dimensional BEC. 

Perpendicular BECs: Another interesting possibil- 
ity is a binary system where the cigar-shaped BEC is 
perpendicular to the plane of the two-dimensional BEC. 
This kind of configuration can be obtained by consider- 
ing the following trapping potential acting on the second 
BEC 

U 2 {r) = W 2 {x,y) + \m 2 u 2 z z 2 , (23) 



Equations (|25j) and ((26)) arc simplified if wc impose cylin- 
drical symmetry in the harmonic potential acting on the 
first BEC, i.e. u)\ x — lo\ v = u>i p , and also in the generic 
potential acting the second BEC, i.e. W 2 (x, y) = W 2 {p) 
with p = (x 2 + i/ 2 ) 1 / 2 . In this way we can choose 
h 2 (x,y,t) = h 2 (p,t) and find the following effective 
BNLSE for fi{z,t) and h 2 (p,t): 
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These are two coupled equations, where the first equation 
is integro-differential in 1+1 dimensions and the second 
is integro-differential in 2+1 dimensions (but in 1+1 vari- 
ables) . 



where W 2 (x, y) is a generic potential in (x, y) plane, while 
the harmonic potential of angular frequency uj 2z along 
the z axis is so strong that the second BEC is a quasi-2D 
in the (x, y) plane. For the first BEC we still use the 
potential of Eq. ([5]) and the variational wave function of 
Eq. (J6j) , but for the second BEC we adopt the variational 
wave function 



V> 2 (r,i) = h 2 (x,y,t) 



e -z 2 /(2ai z ) 
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where a 2z = y hj {m 2 uj 2z ) is the characteristic harmonic 
length of the potential acting on the second BEC along z 
direction. From Eqs. (fTTj) and (fT2]) . using Eq. (|24|) . it is 
straightforward to derive the effective BNLSE for fi(z, t) 
and h 2 (x,y,t) @: 
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C. 3D-2D binary BEC 



Let us suppose now that the first BEC is fully 3D and 
it is under the action of the generic potential t/i(r), while 
the second BEC is trapped by the potential (|23p . Again 
we assume that the harmonic potential of angular fre- 
quency u> 2z along the z axis is so strong that the second 
BEC is a quasi-2D in the (x,y) plane. The first BEC 
is described by the generic wave function i/>i(r,i) while 
the second BEC is modeled by the variational wave func- 
tion of Eq. (f2~4"|) . Following the same procedure above, 
the wave functions ipi(r,t) and h 2 (x,y,t) are found to 
satisfy the following effective BNLSE 
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As in the previous section, these equations are simpli- 
fied if the confining potentials have cylindrical symmetry: 
Ui(r) = U\{p,z) and W 2 (x,y) = W 2 (p). In this way we 
can choose ipi{r, t) = ipi(p, z, t) and h 2 {x, y, t) = h 2 {p, t) 
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and the effective BNLSE become 
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These are two coupled equations, where the first is dif- 
ferential in 3+1 dimensions (but in 2+1 variables) and 
the second is integro-differential in 2+1 dimensions (but 
in 1+1 variables). 



III. DIMENSIONAL REDUCTION BEYOND 
THE EFFECTIVE BNLSE: THE BNPSE 

In Sec. [II] we considered dimensional reduction of the 
GPE for binary BEC in the form of an effective BNLSE 
in reduced dimensions. Such equations are simpler than 
the 3D binary GPE and they could be good approxima- 
tion to the original equations. It is possible to improve 
on the effective BNLSE presented in Sec. [TT| through 
a binary nonpolynomial Schrodinger equation (BNPSE). 
This improvement is obtained by considering a more flex- 
ible variational wave function in place of Eqs. ([5]), (fT5|) . 
and (|24p . where in the transverse direction the wave func- 
tion is taken as the ground state in the respective har- 
monic potential(s) with Gaussian form. In the following 
we maintain the Gaussian form but with a flexible width 
to obtain the required BNPSEs in two cases: (i) The ax- 
ially symmetric 1D-3D case described by Eqs. ([TB"]) and 
(|14[) and (ii) The axially symmetric perpendicular 1D- 
2D case described by Eqs. ([28]) and The other 

cases presented in Sec. [II] can be treated in a similar 
and straight-forward fashion and the respective BNPSE 
obtained. 



A. 1D-3D with cylindrical symmetry 

Here we would like to introduce the next-order correc- 
tion to Eqs. ([TB"]) and (jT4"]) . In this case the confining 
potential ([5]) of the first BEC can be written as 
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We can then use the variational wave function 
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in place of Eq. (O, for the first BEC, where a = a(z, t) is 
its width, which can be quite different from the character- 
istic harmonic length a p ~ yjh/{miuji p ) of the transverse 



confinement. We assume, consistent with the philosophy 
of the reduction scheme, that in Eq. ([33]) the z and t 
derivatives act only on /i . Then the action of the system 
becomes 
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(38) 
(39) 



where we have neglected, in the spirit of the reduction 
scheme, the spatial derivative of a. 

We now extrcmize the action A with respect to f*(z, t), 
^i 2 (p,z,t) and cr(z,t), finding the equations for fi(z,t), 
tp 2 (p, z, t) and cr(z 7 1): 
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Equations gO]) and (04]) are the BNPSE, where the first is 
integro-differential in 1+1 dimensions, the second is dif- 
ferential in 3+1 dimensions (but in 2+1 variables). Equa- 
tion (|42[) is integro-algebraic and determines the width a 
as a function of |/i| 2 and IV^I 2 - If we take a — ai p , the 
harmonic length in the transverse p direction for the first 
component, in place of Eq. (|4"2"]) . we get the previously 
obtained effective 1D-3D Eqs. ([13]) and ([14]). If we set 
the coupling term g 12 = in Eq. (|42p wc obtain the 
model previously obtained in Ref. [B[ for a cigar-shaped 
BEC. 
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B. Perpendicular 1D-2D with cylindrical symmetry 

Here we would like to introduce the next-order correc- 
tion to Eqs. ([2"B"} and ([2"9"} . In this case the confining 
potential is given by Eq. and the variational wave 

function of the first BEC will be taken as (f35|) . The po- 
tential of the second BEC will be taken as 



(43) 



(44) 



U 2 (r) = W 2 (p) + -m 2 ul z z 2 , 
and its variational wave function as 

e -z 2 /(2r, 2 (p,t)) 

where the width r/ = i](p,t) can be quite different from 
the characteristic harmonic length along the z direction 
and the p and t derivatives are supposed to act only on 
h 2 . 

The action of the system can now be written as 
A = J ' L x dzdt + 2tt J L 2 pdp dt + J C 12 d 3 r dt , (45) 

where L\ is given by Eq. (|37|) and 
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i / h 2 



m 2 u} 2 z r] 2 



2 \2m 2 r) 2 



7T 3 / 4 7y(T 2 



\h 2 \ 2 



!J2 



^ 2 | 4 ,(46) 



f 12 \ /.|2|^.|2 



h 2 \ 2 . (47) 



Again we have neglected the space derivatives of the 
widths a and r\. We now extremize the action A with 
respect to f£ , h 2 , a and r\ to get the desired equations 




l/il 2 



e-P 2 /° 2 \h 2 (p,t)\ 2 pdp 



fi 



(48) 



W 2 (p) 



92 



e- z2 l^\h(z,t)\ 2 dz 



(49) 



2 4 

miw lp ~ 



£2 „ r<x> i 

°' = — + f |/i| 2 +4. 9l27 r 1/4 / dpp- 
mi 27r Jo r/ 

-^e-"^^^-^) . (50) 



fr 2 
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m 2 uj 2z ri = 



m 2 2^/27^ 



r 3/4 



a 2 



x e 
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^-2 



(51) 
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FIG. 1: (Color online) Plot of (a) linear and (b) radial densi- 
ties of a binary perpendicular cigar-disk mixture as obtained 
from the 3D GPE Q and Q (3D) and the reduced 1D- 
2D BNLSE ([28} and ((29j) (NL) in dimensionless units for 
Ngi — and Ng2 = 40-7T for different Ng\2- Parameters 
used: mi = 7712 = H = 1, cji p = uj2z = 10, ojiz = ui2 P = 1, 
aiz = «2p = 1, and ai p = a22 = l/\/T0. Both densities are 
normalized to unity. 



Equations fl48J) and 09} constitute the BNPSE, where 
the first one is integro-differential in 1+1 dimensions, the 
second is integro-differential in 2+1 dimensions (but in 
1+1 variables). Equations ([50)1 and (|51[) intcgro-algcbraic 
and determines the widths a and 77 as functions of |/i| 2 
and \h 2 \ 2 . If we take a = a\ p and 7/ = a 2z , the harmonic 
lengths in the transverse directions for the first and the 
second components, respectively, in place of Eqs. ([50} 
and ([Bl"}, we get back Eqs. ([28} and ([29}. If we set the 
coupling term gi 2 = in Eqs. ([50} and ([ST}. Eqs. ([4Jj} 
and (|49p reduce to previously studied models in Ref. [5j 
for a cigar- and disk-shaped BECs, respectively. 



IV. NUMERICAL RESULTS 

We shall present results for the coupled ID- 2D case 
when the cigar-shaped BEC is perpendicular to the disk- 
shaped BEC. The potentials V\(z) and W 2 {p) in Eqs. 
and ([2T)| are taken as 



Vi(z) 



1 



2 2 
■miu) lz z ; 



W 2 (p) 



1 



m 2 w 2p p 2 . 



(52) 



The cigar-shaped BEC, labeled 1, lies along the z axis 
and has trapping frequencies satisfying io\ p jiO\ z = 10. 
The disk-shaped BEC, labeled 2, lies in the (x,y) plane 
and has trapping frequencies satisfying u> 2z /ui2 P = 10. 

To keep the algebra under control, we work in dimen- 
sionless units and set mi = m 2 = h = 1 and take N% = 
N 2 = N. For the cigar-shaped BEC we take oji p = 10 
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FIG. 2: (Color online) The same as in Fig. [TJfor Ngi = 4% 
and Ng2 = 40n. 



FIG. 3: (Color online) The same as in Fig. Q]for Ngi = 4tv 
and Ng 2 = 4007T. 



and uji z = 1; for the disk-shaped BEC we take uj 2z = 10 
and uj 2p = 1- For the cigar-shaped BEC 1, the har- 
monic oscillator lengths are a\ 

V 



y/h/(mu)ig) = 1 and 
For the disk-shaped BEC 2, 



dip = yjh/{rrwj\ p ) = 1/VlO 
the harmonic oscillator lengths are a 2 
l/\/T0 and a 2p = y/K] (mu} 2p ) = 1. 
densities are calculated from the solution of Eqs. ([3]) and 
(® via 



y/h/(muj 2z ) = 
The linear and radial 



|/ x (*)| 2 = 2tt / P dp\Mr)\ 2 , 
Jo 



dz\Mr)\ 2 , 



(53) 
(54) 



to be compared with the solutions of Eqs. (|28[) and Q29J) . 

To solve the coupled GPEs in three and reduced di- 
mensions we employ the split-step Crank-Nicolson dis- 
cretization technique [20| with real and imaginary time 
propagation. We essentially use the Fortran programs 
provided in Ref. [2(| ■ The convergence of the numerical 
result was tested by varying the space (typically ~ 0.05) 
and time (typically ~ 0.001) steps and the total number 
of spatial and temporal discretization points. 



A. Stationary States 

In Figs. Q] (a) and (b) we plot the linear and radial 
densities |/i(z)| 2 and h 2 (p)\ 2 , respectively, obtained from 
the 3D equations (|3) and (gj as well as the BNLSE (gSJ) 
and ([29]) for N g 1 = and Ng 2 = 407r and for different 
Ngn- In this case in the absence of coupling (g\ 2 = 0) 
the two linear densities are identical and the two radial 
densities very close to each other. The BNLSE (|28D and 
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FIG. 4: (Color online) The same as in Fig. [3] obtained using 
the reduced BNPSEs (H) - (E0 (NP). 



(|29|) provide a faithful representation of the densities even 
in the presence of coupling as can be seen from Fig. [1] 

In Figs. [U(a) and (b) we plot the linear and radial den- 
sities |/i(z)| 2 and h 2 (p)\ 2 , respectively, as in Fig. [1] (a) 
and (b), for Ngi = An and Ng 2 = 40ir. In this case even 
for gi 2 = 0, there is some difference between the results 
of the binary 3D CPE and the BNLSE. This difference 
is consistent with the previous studies [H, H| . The agree- 
ment between the densities obtained from the binary 3D 
GPE and the BNLSE continues to be good. The small 
discrepancy between the results of the binary 3D GPE 
and the BNLSE continues as a nonzero 1712 is introduced. 
This discrepancy does not increase as 1712 is increased to 
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FIG. 5: (Color online) (a) Chemical potential of the two 
components /ii and \i 2 calculated using the 3D CPE ([3]) 
and © (3D), the reduced BNPSEs ijlg jl - (f5T |) (NP), and 
the reduced BNLSE (28]) and ([29]) (NL) versus g 2 N/-K for 
ATt?! = 407r,iV 3 12 = 128tt. (b) Chemical potential of the 
two components /ii and fi 2 calculated using the 3D CPE ([3]) 
and @ (3D), the reduced BNPSEs gS]) - (gT]) (NP), and 
the reduced BNLSE (28]) and (29]) (NL) versus giN/it for 
iVg 2 = 400tt, Ng 12 = 128tt. 



moderate values as we see from Fig. [5] In Figs. [3J (a) and 
(b) we plot the linear and radial densities for Ngi = 40tt 
and Ng 2 = 4007T. In this case there is a larger difference 
between the results of 3D and reduced 1D-2D models. 
What is interesting is that the agreement between the 
results of the two models does not deteriorate with the 
increase of gi 2 up to 1287T. 

From the results displayed in Figs. []] — [3J we find that 
the BNLSE (|2"8"|) and ([2"9"]) provide an excellent account of 
the actual state of affairs, when compared with the full 
3D binary GPE (3J and (3]), specially for small nonlin- 
earities Ngi and Ng 2 . 

We next attempt an improvement over the BNLSE 
([28]) and (29]) of Sec. [TT]by applying BNPSE of Sec. [En] 
In other words, we attempt a solution of Eqs. (T4"8"]) — ([5T]) , 
in place of Eqs. (f2"B"j) and (f2"9"]) . to reduce the discrepancy 
noted in Figs. [3] (a) and (b) between the results of the 
3D and BNLSE. The results so obtained are plotted in 
Figs. H (a) and (b) for Ngi = 40tt and Ng 2 = 400tt, 
which show a significant improvement over the results in 
Figs. [3J (a) and (b). The large discrepancy between the 
results of the 3D binary GPE and the BNLSE for the 
linear density of component 1 has virtually disappeared 
in Fig. @] (a) and that for the radial density of compo- 
nent 2 is significantly reduced in Fig. |H(b). This shows 
that although the results from the effective BNLSE of 
Sec. [II] arc quite good for small nonlincarities, significant 
improvement can be obtained from the BNPSE of Sec. 
IIII1 especially for larger values of nonlincarities. 
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FIG. 6: (Color online) Evolution of axial and radial rms 
sizes of the cigar- and disk-shaped components of the bi- 
nary BEC after a sudden reduction of the radial angular fre- 
quency of the disc-shaped component by 5%, for the non- 
linearities (a) Ngi = 4ir,Ng2 = 4071,^312 = 16n and (b) 
Ngi = 407T, Ng 2 = 4007T, Ng 12 = 128n. The cigar-shaped 
BEC is perpendicular to the disk-shaped BEC. 



Finally, we consider the chemical potential fi\ and [ii 
of the two components calculated using the 3D binary 
GPE (3J and ©, the BNPSE (|3g]l - (|5T|) as well as 
the BNLSE (HI) - (EH)- The equations for the chemical 
potentials are obtained by replacing the time derivative 
ihdt by /ii and /i2, respectively, in the dynamical equa- 
tions for component 1 and 2, respectively. The chemical 
potentials are plotted in Fig. [5] (a) versus giN/ir for 
Ng 2 = 400tt and Ng 12 = 128tt. In Fig. [5] (b) we plot 
the chemical potentials versus g 2 N/^ for Ngi = AOtt and 
Ngi 2 = 128-7T. The chemical potentials obtained from the 
BNLSE (5HJ) - (23) are good but those calculated from 
the BNPSE (48]) - (51]) are better approximations to the 
results obtained from the 3D binary GPE (3J and (5J) in 
all cases. 



B. Dynamics 

Now we investigate if the reduced equation presented 
here can satisfactorily describe the dynamics of the per- 
pendicular 1D-2D system. For this purpose first we con- 
sider the coupled system presented in Figs. [2] (a) and 
(b) with Ngi = 4tt, Ng 2 = 40tt, and Ng 12 = 16tt. After 
the stationary condensate is created during time evolu- 
tion of the numerical routine using real-time propaga- 
tion, we reduce the radial angular frequency of the disk- 
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shaped quasi-2D BEC suddenly by 5% (the new angular 
frequencies being uj\ p = cu 2z = 10, u)\ z = l,tu 2p = 0.95,) 
and study the subsequent dynamical evolution of the sys- 
tem. The relevant dynamics is best illustrated by plot- 
ting the evolution of the rms axial size of the cigar-shaped 
quasi-lD BEC and the rms radial size of the disk-shaped 
quasi-2D BEC in Fig. [5] (a) and (b) for two sets of non- 
linearities exhibited in Figs. [2] and [3] Here we show 
the results of the 3D GPE © and g]) (3D), the reduced 
BNPSEs gHl) - (ED) (NP), and the reduced BNLSE (gSJ) 
and (f2"9"|) (NL). For both sets of nonlinearities illustrated 
in Figs. ED (a) Ng x = in,Ng 2 = 40tt, Ng X2 = 16tt and 
(b) Ngi = 40tt, Ng 2 = 400tt, Ng l2 = 128tt, the NP re- 
duction produced far better results compared to the NL 
reduction. In both cases the angular frequency of the 
quasi-2D disk-shaped BEC was modified and both the 
NL and NP schemes as well as the full 3D calculation 
exhibited sinusoidal oscillation for the axial size of the 
disk-shaped BEC right from t = 0. It took some five 
units of time for similar sinusoidal oscillation to initiate 
in the unperturbed cigar-shaped component due to the 
non-zero coupling term g\ 2 . The oscillations in rms sizes 
in the NP scheme were always in phase with those of the 
full 3D calculation. But for larger nonlinearity and at 
large times the oscillations in rms sizes in the NL scheme 
are not always in phase with those of the full 3D calcu- 
lation. The differences between the rms sizes of the NL 
and NP schemes compared to the full 3D calculation were 
consequences of the same differences, existing at t = 0, 
in the corresponding stationary results in the absence of 
any perturbation in the angular frequency. 

V. SUMMARY 

In this paper we derived and studied simple binary re- 
duced non-polynomial Schrodingcr equations (BNPSEs) 
for the description of a binary BEC where the two compo- 
nents subject to distinct trap symmetries belong to two 
different spatial dimensions. We considered three possi- 
bilities: the first component is in 3D and the second in (a) 
ID (cigar shape) or (b) 2D (disk shape), and (c) the first 



component is in ID and the second in 2D. The ID and 
2D configurations were achieved in an axially-symmetric 
setting with a strong harmonic trap in the transverse ra- 
dial and axial directions, respectively. The BNPSEs were 
obtained by a Lagrangian variational scheme after ap- 
proximating the spatial wave-function component in the 
transverse direction by a generic Gaussian function. A 
simpler set of reduced equations, the BNLSE with power- 
law nonlinearity, were obtained with a simple ansatz for 
the 3D wave function where the wave-function compo- 
nent in the transverse direction is taken to be the ground 
state in the corresponding harmonic trap. 

We tested the accuracy of the reduction scheme above 
by numerically solving the 3D binary GPE in the coupled 
1D-2D case where the cigar-shaped component is in the z 
direction and the disk-shaped component lie in the (x, y) 
plane. We calculated the density profiles and the chem- 
ical potentials of the two components. For smaller val- 
ues of inter- and intra-species nonlinearities, the BNLSE 
([2"5]) and (|2"9"|) with power-law nonlinearity produced re- 
sults for density and chemical potential in good agree- 
ment with the binary 3D GPE. For larger nonlinearities 
the BNPSE ([48]) - ([51]) produced significant improve- 
ment over the BNLSE J2SD and flU]). We also studied 
small oscillation of the binary 1D-2D BEC, initiated by 
a sudden change of the radial angular frequency of the 2D 
component. For larger nonlinearities and large times, the 
BNPSE (gHl - dH]) produced better results for dynamics 
compared with the BNLSE (g5]) and (JUJ). To conclude, 
the BNLSE ([28]) and ([29]) as well as the BNPSE ([48]) - 
(|51[) are shown to be very useful to investigate a binary 
1D-2D BEC, where the cigar-shaped BEC is perpendic- 
ular to the disk-shaped one. Similar reduced equations 
are derived for other symmetries, such as the 1D-3D and 
2D-3D cases as well as the 1D-2D case where the cigar- 
shaped BEC is in the plane of the disk-shaped one. 
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